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Abstract 

The renormalization group method has been adapted to the analysis of the long-time 
', behavior of non-linear partial differential equation and has demonstrated its power 

' in the study of critical phenomena of gravitational collapse. In the present work we 

apply the renormalization group to the Einstein equation in cosmology and carry 
out detailed analysis of renormalization group flow in the vicinity of the scale invari- 
ant fixed point in the spherically symmetric and inhomogeneous dust filled universe 
model. 

> ' 

O 

Qs, • I. INTRODUCTION 

O 



PACS numbers: 64.60.Ak, 98.80.Hw 



Recently the renormalization group (RG) idea is applied to study the long-time asymptotics of the non- 
hnear partial differential equations . The RG transformation there is integration of the equation up 
to a finite time followed by a rescaUng of the dependent and independent variables. The RG transforma- 
tion together with the original differential equation gives a RG equation. Using the RG transformation, 



5-H . the problem on an infinite time is reduced to the problem on a finite time. A fixed point of the RG 

transformation corresponds to a scale invariant solution of the differential equation. We can obtain the 
long-time behavior of the equation by studying the flow around fixed points. 

As an application of this RG method to the system of gravity, Koike, Hara and Adachi |^ analyzed 
^ , the Einstein equation to understand the problem of critical behavior of black hole mass in gravitational 

' collapse found by numerical study [Q. A pedagogical exposition of the RG method in the deterministic 

system is given by Tasaki in a simple but very illustrative example of motion of a point particle in the 
Newtonian gravity. 

Here we apply the RG to the Einstein equations in the cosmological situation. For simplicity, we shall 
consider only two cases. One is a homogeneous and isotropic universe filled with a perfect fluid and the 
other is a spherically symmetric universe filled with dust. We shall study the flow near the fixed points 
of the RG equations, which have self-similarity. 

The astronomical observations indicate that the present universe has the hierarchical structure such 
as galaxies, clusters of galaxies, and super clusters, and that the two point correlation function of the 
galaxies and of the clusters of galaxies can be expressed roughly by a single power law The scale 
invariant Harrison-Zel'dovichi spectrum for the primordial density perturbation has been successful in 
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the study of structure formation of the universe |7|] . These suggest that the present universe has some 
self-similarity and that the scale invariant solution plays an important role in cosmology. 

In this paper we apply the renormalization group method to the Einstein equations in the cosmological 
context. In Sec. II, we illustrate the application of the RG method to the heat equation with the nonlinear 
term. In Sec. Ill, we apply the RG method to the Einstein equations. Section IV is devoted to the summary 
and discussions. 



II. RENORMALIZATION GROUP TRANSFORMATION HEAT EQUATION WITH 

NONLINEAR TERM 



In this section, we review the RG method for nonlinear partial differential equations 
consider the heat equation with the nonlinear term as a simple example: 



du{x, t) 
dt 



1 



i"{x,t) + \u'^{x,t)] , 



First we 



(2.1) 



where the prime denotes the spatial derivative and A is a coupling constant. The equation (2.1) has scale 
invariance under the following scale transformation: 



X 

t 

u{x, t) 



Lx, 

L^u{Lx,L'^t), 



(2.2) 



where L is a par amet er of the scale transformation and is taken to be larger than 1. Namely, if u{x, t) is 
a solution of Eq.(|2.lD, the scaled function. 



u {x,t) = L^u{Lx,L^t), 



(2.3) 



is also a solution of Eq.( |2.1| ). We can thus obtain a one-parameter family of solutions, provided that 
u(x, t) is a solution. 

Here we define the RG transformation 7?.^ of a function of x by 



TZlu{x, 1) ~ u (x, 1) 



(2.4) 



In short, the TZj^ is a map from a set of initial data to another. It is convenient to take the initial time 
to be t = 1. The RG transformation TZ^ has a semi-group property: 



7^r,^-l oUr 



(2.5) 



Letting t — 1 and then — t \n Eqs.(^.3|) and (2.4), we can express an arbitrary solution u{x,t) as an 
initial data u(-, 1) transformed by RG transformations: 



u{x,t)=t-^ ^\ ^ (xt-i/2,1). 



(2.6) 



The large L means the late time. Repeating the RG transformation ( p.4[ ), we can see the long-time 
behavior of the solution u{x,t) in Eq.( ^.lD . 
Denoting L — e'^, we have from Eq.(p.3D that 



du ^du JL) du 
— ^ L— ^2u +xu' +2 — . 
dr dL at 



(2.7) 



Using the original partial differential equation (2.1) we have 



d u (L) 

— 2u 

dr 



(L) 



xu' + u" . 



(2., 
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This is the equation satisfied by the scaled function u . which we caU the RG equation. We note that the 
equation ( |2.8| ) has no expHcit scale L dependence because of the scale invariance of the original equation 

We investigate the fixed point of the RG equation ( |2.8| ). The fixed point u* is defined by 



TZlu* = u* 



(2.9) 



for any L > \. This condition means that the field profile is unchanged after time evolution followed by 
suitable rescaling. In general this condition is equivalent to 



du* 
1h 



= 0. 



From Eq.( p.8D , u* satisfies the following equation, 

2u* + \u''^ + xu" + u*" ^Q. 
In the homogeneous case, we can easily obtain the fixed points: 

2 



= and — 



A' 



(2.10) 



(2.11) 



(2.12) 



To investigate the character of the fixed point Eq.(2.12), we consider the linear perturbation around 



the fixed point Eq.(2.12). The perturbed quantity 5u is defined by 

(L) W 

u — u* + Su, 



(2.13) 



where Su is assumed small. Substituting Eq.(2.12) into Eq.( |2.8| ) and neglecting the second order term 

(L) (L) 

SvF' , we obtain the linearized equation for 5u: 



dSu 



We require the boundary condition 



(L) (L) (L) 

2du + xSu' + 5u" {u* = 0) 

(i) (i) (L) 

-2Su + xSu' + 6u" (u*=-f) 



(i) 
Su 



(2.14) 



{\x\ — *oo). 
We are going to find the normal modes with the ansatz: 



Su — f{x)e 



(2.15) 



(2.16) 



where f{x) is a function to be determined below and w is a constant. From Eqs.(2.14) and (2.16), we 
have 



/" - xf - (CO - 1)/ = {u* = 0), 



f"-xf-{u; + 3)f^0 {u* 



The regularity at x = and the boundary condition at = oo imply 

fix) = H^ix), 

1-n {u* = 0) 
'""^ -3-n {u* = -j:) ' 



(2.17) 

(2.18) 
(2.19) 
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where Hn(x ) is th e Hermite polynomial and n = 0, 1, 2, ■ 



From Eg. ( 2. IE ), u* = — 2/A is an attractor because all w's are negative. On the other hand, u* = 



has only one relevant mode {n — 0). 



We can discuss the long-time behavior of a solution of th e non linear diffusion equation (2J), if u{x, 1) 



is sufficiently close to the self-similar profile, u* . From Eq.(2.12), we obtain the two self-similar profiles. 
Suppose the initial spatial profile of Su is expressed as a superposition of the normal modes i/„e~^ 



As we have seen from Eq.(2.19), if u(x, 1) is sufficiently close to the fixed point u* ~ —2/A, the solution 
approaches to —2t~^/X in course of time because all modes of perturbation are irrelevant. On the other 
hand, there is the only one growing mode ( n = ) of perturbation around the fixed point u* = 0. As 
time goes on, the behavior of the solution near this fixed point is dominated by the relevant mode n — 0. 
This relevant mode corresponds a gaussian distribution t~^/^e~^ /(2t) 

From this instructive example, we see that if the perturbation around the fixed point has a finite number 
of relevant modes or no relevant modes, we have some prediction power for the long-time behavior of 
nonlinear partial differential equation. 

III. RENORMALIZATION GROUP FOR EINSTEIN EQUATION 

In this section, we apply the RG method, which is explained in the previous section, to the Einstein 
equations. 

Take a synchronous reference frame where the line element is 

= -df +j,j{x'',t)dx'dx^. (3.1) 

Throughout this paper Latin letters will denote spatial indices and Greek letters spacetime indices. The 
matter is taken to be a perfect fluid characterized by the energy-momentum tensor 

Tf,i, = (p + p)u^u^ + pg^^, (3.2) 

where p, p, and are pressure, energy density, and four-velocity, respectively. We assume that the 
equation of state of the fluid is 

p={T-l)p, (3.3) 
where F is a constant. The Einstein equations are 

%^2K,,, (3.4) 

k,, = -^R,, ~ KK,, + 2K\Ki, + ^ [2Tu,Uj + (2 - F)7y] , (3.5) 

2(1 -I- Fu;m') 

KTpu, = , ^ {Kr~K^^, (3.7) 



uiw 



where Kij is the extrinsic curvature, '^Rij is the Ricci tensor associated with jij, and k = 8ttG. A dot 
denotes the derivative with respect to t, a semicolon denotes the covariant derivative with respect to -fij. 

Hereafter we consider the RG transformation for the dynamical variables ^ij and Kij . In the following 
subsections, we investigate the two cases; the one is a homogeneous and isotropic universe, the other is 
a spherically symmetric inhomogencous dust universe. 



A. Homogeneous and Isotropic case 

We consider the homogeneous and isotropic universe as a simple case. This case is rather trivial because 
the fleld equation becomes an ordinary differential equation. Nonetheless, this gives a nice warming up 
model to familiarize us to the RG approach to the universe. In this case, the spatial metric is written by 



4 



1 



k(t)r 



2 ■ 



(3.8) 



where = S ^x'^ x^, and a(t) and k{t) are the functions of time t to be studied. Substituting Eq.(3. 
into Eqs.(O) - we get 



a = aH, 
H = -3i?2 
fc = 0, 



2k 

2 



-Kp, 



Kp — 3 



where H is the Hubble parameter. From the conservation law of the energy density, we have 



-3r 



Kp — Ala 

where M is an arbitrary constant. 

First, we consider the following scale transformation. 



(3.9) 
(3.10) 
(3.11) 
(3.12) 

(3.13) 



a{t) 

m 



t — > Lt, 



a\t) = L-'/^'^^aiU), 



(i). 



k (t) = L 



^ r2(3r-2)/(3r) 



k{Lt), 



(3.14) 
(3.15) 



where L is a parameter of scale transformation and larger than 1. From Eqs.( 3.10 ) and ( 3.13 ), under this 
scale transformation the variables H and p are scaled in the following way, 

(L) 

H{t) = LH{Lt), 



H{t) — > H{t) = LH{Lt), (3.16) 
Pit) ^ ^p\t) = L'piLt). (3.17) 

The equations of motion (3.9)-(33J) are invariant under the scale transformation Eqs.(3.15)-(3.17) and 
equivalently the scaled variables satisfy the original equation. 
Second, we define the RG transformation TZl'- 



7^La(l) = (1), 
Letting t = L, we have formulas 



TIlH{1) = Hil), 7^Lfc(l) = k (1). 



ait) = t2/(3r) (*) 

H{t)=t-'H (l), 
k{t) - t-2(3r-2)/(3r) 

which we shall use later to see the long-time behavior of a, H, and k. 

Third, we derive the RG equation. Letting L = , the infinitesimal transformation of a , , and k 
with respect to r is 



(3.18) 

(3.19) 
(3.20) 
(3.21) 



a a 

(L) 

dH 
dr 

d k 
dr 



2 (i) 
a 

3r 



a a 



H 



dH 

dt ' 



2(3r-2)(i) dk 
k + 



3r 



dt 



(3.22) 
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Using the equations of motion (3.9), (3.10), and (3.11), Eqs.(3.22) can be rewritten as 



^^^^ 
a a 

(L) 

dH 
dr 

(L) 

d k 
dr 



2 {L) (L){L) 

a + a H , 

3r ' 



H 



2-3r (L) 

K P 



2(3r- 2)(i) 
3f ^ ■ 



(3.23) 



These equations (3.23) are the RG equations. 

Here we investigate the fixed point of the RG equations. The fixed point (a*, H* , k*) is defined by 



7^La* = a*, TZlH* 
The above conditions can be rewritten as 



H*, 



(3.24) 



da* 



0, 



(L) 

dH* 



0, 



From Eqs.( 3.12 ) and ( 3.23 ), the fixed point is 



dk* 
~d^ 



(3.25) 



MIT' 



^ -3f' 



0, Kp* = 3H* 



(3.26) 



This fixed point corresponds to a fiat Friedmann universe. 

Note that if F is taken to be 2/3, there is another fixed point where k* is non-zero. For the non-zero 
k* case, the term of the spatial curvature can be absorbed into the term of the energy density of matter 
because the dependence of the scale factor on each term is the same. Thus the k* = case includes the 
non-zero k* case. Hereafter we concentrate on only the fc* = case. 

In order to study the flow in the RG around the fixed point, we consider the perturbation around the 

(L) (L) (L) 

fixed point. The perturbed quantities 5a, 5H, and Sk are defined by 



a = a 



(L) 

5a, 



H =H* 



(L) 

dH, 



(L) (L) 

k =k* + 6k, 



(3.27) 



(L) (L) (L) ^ ^ 

where 5a, 6H, and 5k are assumed small qualities. From Eq.( 3.23| ), the perturbed quantities satisfy the 
linearized equations. 



d5a 

(i) 
d5H 

(L) 

dSk 



a*5H, 



(L) 

^5H 



2a* 



2(3r-2)(4) 

3F 



where we neglect the second order term Sa"^ and 5H^ and use the linearized equation of Eq.( [3.12 ): 



(i) 

K.5p 



4 W 
-5H 



(3.28) 
(3.29) 
(3.30) 

I ; 

(3.31) 



Substituting Eq.( |3.28D into Eq.(p.29|), Sa satisfies 
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We solve Eq.(|3.32|) 



(P5a dSa 3r „ 



(-^') 2(3r-2) 

= /le-" + /2e^^^ (3.33) 



where /i and /2 are arbitrary constants. From the solution ( 3.33| ), we can see the flow in RG around the 



fixed point. If 3r — 2 < 0, this fixed point is an attractor. On the other hand if SF — 2 > 0, there is a 
single relevant mode. Note that in the case 3F — 2 > 0, the matter we consider satisfies the strong energy 
condition. 

From the flow in the RG around the fixed point, we can see the long-time behavior of the homogeneous 
and isotropic universe. If 3r — 2 < and setting the initial profile in the vicinity of the fixed point 
{a* , H* ,k*), the spacetime will approach to the fiat Friedmann universe a {t) = a*t2/3r^ On the other 
hand, if 3r — 2 > 0, the spacetime will deviate from the fiat Friedmann universe because there is a relevant 
mode Sa{t) = f^f<^3r-2)/3r^ 

In the context of the usual cosmological perturbation around a fiat Friedmann universe, /i mode 
corresponds to the decaying mode and /2 mode corresponds to the growing mode, which implies the 
gravitational instability, because the matter should satisfy the strong energy condition. 



by 



B. Spherically symmetric case 

We consider the spherically symmetric inhomogeneous case. In this case, the spatial metric is written 



7y (2:^ i)dxMx^' A^{r, t)dr'^ + B'^{r, t){d0'^ + sin^ Od^''^). 



(3.34) 



Namely, 7rr = [r, t) , 795/ = (r, t) , and 7^^ = sir? 9 while the other components of the spatial 
metric vanish. As a s imp le case, we investigate the universe filled with dust, i.e. F = 1 and we can set 
= in Eqs.(3.4) - (3.7). The Einstein equations are 



1^3 



2K, 



2KaK\ - KK, 



Kp^-l^R + K^-K^Kl] 



Here we consider the following scale transformation. 



(3.35) 
(3.36) 

(3.37) 
(3.38) 



r 
t 

Irrir, t) 

lee{r, t) 
K„{r, t) 
Kee{r, t) 



Lr, 



-2-2q 



7„(Lr, L'^t), 



(r,t) = L-2"7ee(ir, L'^t), 
Krr{r,t) = L^-°'K„{Lr, L°'t), 

(L) 

Kee{r,t) = L~"Kge{Lr,L''t), 
^'p{r,t)^L^"p{Lr,L"t), 



where a is an arbitrary constant because the coordinate transformation, r ~> 
a/f3 for a in Eqs.( 3.39 )-( 3.43 ). Without loss of generality, we take a to be positive. Because of the scale 



(3.39) 
(3.40) 
(3.41) 
(3.42) 
(3.43) 

yields substitution of 
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F — h — K • W W 

invariance of the Eins tei n equ ations Eqs.( p.35 ) - (3.38), the scaled variables jrr, jee, K„, Kqq, and P 

also satisfy Eqs.(^^ - ( p8| ). 



We derive the RG equation. Letting L 
respect to r is 



e"^, the infinitesimal transformation of 7^ and iiTy with 



ajrr 

Ik' 



dr 

(L) 

dT 

(L) 

dKge 
dT 



(L) (L) W 

Z(i — an 

(L) (L) W 

-2a-fee + rdr-fee + ^.aKee, 



(2 - a)Krr 



(L) 



-1 Rlrr 



(L) (L) (i) (L) 
2KrlKl- KKrr 




K ^rr 



W (L)^ 



-aKaa 



rdrKaa 
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(i) 

'Ron 



2KeiK'-KKgg + - K 



(3.44) 
(3.45) 
(3.46) 

(3.47) 



where Rij is th e Ricc i ten sor associated with 7^^. In th e der iva tion o f Eqs.(3.44) - (3.47), the equations 
of motion (3.35) and ( 3.36| ) are used. These equations ( ^.44 ) - ( ^.47 ) are the RG equations. 



W (f) (L) 



(i) 



In terms of A and B {A^ ~ ^„ and — 76(e), the RG equations (3.44)-(3.47) read 



dr^ 



-r^A" 



1 dB 



rB' 



2(«-l) + 77T-^-T7T 



(i) dT 

B 



(2a- 3)rA' 



(L) 

rB' (L) 

~^ 

B 



B 



dA 
dT 



rB' 



2r- 



dA' 



{L)iL) 

(L) r^A'B' 

A - 



(L) {L)(L){L) (L)iL) 

a^{A^-2A' BB' - AB'^ 



(L)(L)(L) 

2 A BB") 



d^ B 
dr^ 



-r'^B" 



2a 



iL)(L) 

2A2B2 
J_dB_ _ 

(L) dT 

2B 



52 



B 

(L) (L) 

A dB 

2 dT 



4a + 2(i) 
-7: A 



{L)(L) (L)(L) (L),L) 

A B -r AB' + rA' B 



rB' 
B 



(i) 
dB 



2r- 



dB' 



{2a - l)rB' 



(L) 

2B 



(L) (L) 

a^jB'^^A^) 
2A2B 



(3.49) 



where the prime denotes the derivative with respect to r. From Eqs.( |3.37 ) and (3.38), we obtain 
2 



K P 



3a - 



q;2 a 
2 



1 dB 



(L) dT 

B 



(i) 

rB' 
B 



dA 



dT 



a2 B 



2(2a- 1) 



(L) 

^rA' 
A 



(L) 

^rB' 
B 



(L) 

1 dB 



(L) dT 

B 



dB 



dT 



A^B^ 



(L) 

2 rA' 

a~ 
A 



(L) (L)(L) (L) 

2{2a-l)rB' 2 r'^ A' B' 1 r'^B''^ 



B 



a2 (L)(L) 

A B 



52 



A^ 



2A' B B' - AB'^ -2A B B" 



(3.50) 
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B'^- A- rA' B' +r AB" = Q. 

dr dr 



(3.51) 



Letting t = L°', the original variables A{r,t) and B{r,t) are expressed by the scaled variables A and 



B: 



(t 



B{r,t) ^t B' 



Here we investigate the fixed point of the RG equations (3.48) and (3.49) defined by 



dA* 
'dT 



= 0, 



dB* 



= 0. 



At the fixed point, 



,l/o 



A{r,t) = A* (rt-i/", 1) = x ( function of rf-^/" only ) 



and 



B(r, t) = t B* (rt-i/", 1) = t x ( function of rt"!/" only ) 



(3.52) 
(3.53) 

(3.54) 

(3.55) 
(3.56) 



is a self-similar solution. 

In the spherically symmetric spacetime filled with dust, the general solution of the Einstein equations 
is the Tolman-Bondi solution (Al) - (A2) in the Appendix. Therefore we can obtain the fi xed point from 
the Tolman-Bondi solution with self-similarity rather than solving the equations (3.54) directly. The 
precise form of these are 



For c = Q 



ar°/3-i(l-3pr° 

A (fA) — 1 — 

3(l-]5r«)3 

B*{r,l) =r"/3(l-pr")i, 
4 



3(1 -pr")(l - 3pr")' 



For c > 



For c < 



A*ir,l) 
B*{r,l) 
sinh rj — rj 
Kp*{r,l) 

A*(r,l) = 



(1 + c)i/2r 
2 

— r"(cosh7y — 1) 
9c 

0^3/2 



— (coshTy — l)r' 
9c 



« _ ^1/2. 



sinh r/ 
cosh T] ~ 1 



9c' 



'(cosh^-l)2 ^(cosh,y-l)-cV2_|Bk^ 



sm77 



(1_ |c|)l/2j. 

i?*(r,l) = -^r"(l-cosr]), 
9|c| 

9|c|3/2 _ 
ri-smr]= — - — [r -p), 



up* (r, 1) 



— (l-cos77)r"-|c|i/2 
9|c| 1 — cosry 



9c2 



r"(l-cos,7)2 2^(l-cosry)-|c|i/2_^ 



3.57) 

3.58) 
3.59) 

3.60) 
3.61) 
3.62) 
3.63) 

3.64) 
3.65) 

3.66) 
3.67) 
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where c and p are constants. 

The constant c can be interpreted as the total energy of the universe in the analogy of the Newtonian 
mechanics. By the signature of the constant c, these fixed points are classified into the following three. 
The universe with c = is similar to the fiat Friedmann universe. The universes with a positive c 
or a negative c are similar to the open and closed Friedmann universes, respectively. Especially when 
c = p = 0, the above fixed point coincides with the fiat Friedmann universe and the spacetime becomes 
homogeneous. 

In the context of the RG, we can treat the time evolution of the field variables as the map from a 
set of initial data to another. If the initial data is taken to be the above fixed point A*, B* , and p* , 
the spacetime will evolve into the Tolman-Bondi solution with self-similarity. In the cases of c = and 
c > 0, the fixed point is not regular if p > 0. For c < 0, the fixed point has singularities irrespective of 
the signature of p because the spacetime is similar to the closed Friedmann universe and will recoUapse. 
Since we should set a regular initial data in the physical situation, we investigate only the case of c = 
(p < 0) and c > (p < 0) where the fixed point is everywhere regular. Note we exclude the case of 
c = p = because we have already studied it in the previous subsection. 

To study the behavior of the flow in the RG around the fixed point, we consider the linear perturbation 
around the fixed point. The linear perturbation around the self-similar Tolman-Bondi solution has also 
been discussed by Tomita in a different context from ours Jl3t . For simplicity, we concentrate on the 

(L) (L) 

spherical modes of linear perturbation. The perturbed quantities 5 A and SB are defined by 



7„ = A*2 + 2A*5A, 

(L) „ (i) 

7ee = 2B*SB. 



(3.68) 



We assume the spatial metric variables for S A and S B the following form. 



(i) 

SA = a(r)e'' 
SB = h{r)e'^ 



(3.69) 



[L) 

The perturbed quantities a(r) and 5 P are expressed by b{r)\ 



a 

A* ^ B' 
S P 



2(1 + c) 
9{uj + a) 



u b b' 
4a '^B^^B^ 



(3.70) 
(3.71) 



where Ci and C2 are arbitrary constants (see Appendix). 

As for the spherical modes of the perturbation, we can easily obtain the solutions. 
For c = 0: 



6(r) = 



(3.72) 



where C3 is another arbitrary constant. The density contrast is 
9 



S P 



20a 



r-2a/3(i _ pr")2/3(l - 3pr")-i [3t^ -f a - 3(tJ -f 3a)pr"] 



o o 

-— C2pr"(l - 3pr")"i + -C3r-"(1 - pr")-i(l - 3pr")"i [w 2a - (w + 3a)pr"] 
2a a 



(3.73) 



In the expression for the linear perturbation Eq.(3.73) there are three terms corresponding to ci,C2, and 
C3 so that there should be a gauge mode hidden in Eq.( 3.73| ) because the number of physical modes has 
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t o b e two. Actually there remains a gauge freedom corresponding to the coordinate transformation of r 
([A^). The gauge mode is given by 

p* p* 

= e'^^r'^f [2apr"(2-3pr")(l-pr")-i(l-3pr")"i] , (3.74) 

with / being an arbitrary constant. Because we should fix the freedom of gauge we choose / so that 

{S p + S p g)/p* behave as nicely as possible at r = oo because we are interested in the perturbation 
modes which are finite at r = oo. 

We use the following condition as a convenient gauge condition. 

[L] (L) (L) 

We fix the gauge mode by the above condition and obtain the physical perturbation 6 p = S P + 5 p g 
as follows. 



p* 



e^'^r" {Ai(1 - 3pr")-i p-2/3^-2a/3(^ _ pr")^/^ [3uj + a ^ 3{u} + 3a)pr"] 
+ (tj + 3a)pr"(2-3pr")(l-pr")-i] + A2pr"(l - pr")-i(l - 3pr")-i} , (3.76) 



where 



Ar = -4^c„ (3.77) 
A. = -|(c.-±c3). (3.78) 

Since we consider only the case of p < 0, the coordinate r can be taken from to oo. We demand 
. . (i) 

the regularity condition that S p / p* should be finite at the boundary, r = and r = oo. This condition 
implies that Ai modes with 2a/3 < uj < a and A2 modes with —a < uj < a are allowed. The mode 



with UJ ~ corresponds to change of p in the self-similar solution Eqs.(3.5!;) and p* remains constant 
independent of t in the direction. Although this fixed point is not a repeller, it has many relevant modes, 
Ai with < u! < a and A2 with < w < a. Note that a suitable linear combination of the Ai and A2 
modes will have an asymptotic behavior w j.'^-'^a at r = 00. For such modes, 2a/3 < uj < 2a is allowed. 
These modes which satisfy the regularity condition are not discrete but continuous. We suppose that 
this special feature arises because our matter is assumed to be dust. To summarize, the possible value of 
w ranges from —a to 2a. The flow of RG in the vicinity of the fixed point is shown in Fig.(|^). 

In the case of c = p = w here t he fixed point corresponds to the flat Friedmann universe, the C2 



modes in the density contrast ( 3.73 ) are gauge modes. The regularity condition implies that ci mode 
with UJ = 2a /3 and C3 mode with uj — —a are allowed. This result corresponds to the homogeneous 
and isotropic case in the previous section. Compared with the usual cosmological perturbation in the 
synchronous comoving reference frame, this result appears to be strange because the only spherical modes 
of the linear perturbation allowed are constant in space. However, the time coordinate used in this RG 
method is different from the usual cosmic time coordinate, the solution allowed by the regularity condition 
in each case does not coincide in general. Moreover, in the homogeneous universe, there is no non-trivial 
characteristic profile of field variables. If the fixed point is a homogeneous universe, the RG method may 
have no advantage since the RG approach respects the self-similar profile. But if the fixed point is an 
inhomogeneous universe, we believe that the RG method may be useful. 

Compared with c = 0{p < 0) case, the value of uj allowed in the case {p — 0) is the lower limit in the 
case {p < 0). The effect of nonlinearity of gravity makes the growth rate of the density contrast large. 

For c > 0: 
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6(r) = r'^ 



9c2 



Cl 



2(cosh77 — 1) — 



3 sinh 77(sinh r] — rj) 



cosh 77 — 1 



— C2 
2c 



cosh 77 — 1 — 



sinh77(sinh7] — 77) 
cosh ?7 — 1 



c^/2 sinh 77 
cosh 7y — 1 



C3 



(3.79) 



The density contrast is 



S P 



-(cosh 77 — 1) 



cosh rj— 1 



9c2 



Cl 



(cj + 3a)r" 



2(cosh77 — 1) — 



3 sinh 7/(sinh r] — r]) 



C2 

2c 



+Cl/2C3 



cosh 77 — 1 
sinh 77 (sinh 77 — rj) 
cosh 77 — 1 
^ sinh 77 9q:c2/'^(2 cosh77 + 1) 



2a(cosh?7 — 1) — (uj + 3a) 
{uj + 3a)r' 



27c^/2a 
2 

9c2/3 



sinh?; (2 cosh 77 + 1) (sinh 77 — 77) 



cosh 7/ — 1 



(cosh?/ — ly 



{uj — a) sinh 77 a{2coshri + 1) (sinh 77 — ?;) 



cosh 77 — 1 



(coshry — 1)^ 



cosh 77 — 1 
The gauge mode is given by 

p* p* 



2(cosh77- ly 



(3.80) 



¥(cosh77-l)-^^:^ 



4 , , Ac^/^smhri 9c^r^"(2 coshri + 1) 

— r (cosh 77 — 1) \ 

9c cosh 77 — 1 2(cosh77— 1)^ 



(3.81) 



Similarly to the case of c = 0, we demand the regularity condition at r = and r = 00. 

(i) case 1 {p <0): 

We use the following convenient gauge condition, 



1 ( UJ + 3a 



-Cl 



2-3 



sinh ?7o (sinh 7/0 - t/q) 



:C2 



2a — {uj + 3a) 



(cosh 770 — 1)^ 
sinh ?7o (sinh 7/0 - t/q) 
(cosh 7/0 — 1)2 



9c2/3(cj + 3a)(sinh 770 ) 
(cosh 770 — 1)2 



-C3 



(3.82) 



where rj = rjo corresponds to r — > 00 and 7/0 is thus given by sinh 770 — 770 = —9pc'^^^/2. 

, -1 W (L) (L) 

By using the above condition Eq.( |3.82| ), we obtain the physical perturbation 6 p — 6 P + S P g as 
follows. 



(i) 
5p 



a 



^(coshT?- 1) - £t!llinhv 



{ 



At 



3(w + 3a)r"(cosh7y- 1) 



cosh ?/ — 1 

sinh r/o (sinh 7/0 — i]q) sinh 77 (sinh 77 — 77) 



9c3/2 



4cj + 9a - 6{uj + 3a) 



(cosh 770 — 1)2 
sinh 770 (sinh ?7o - ?7o) 



(a; + 3a)r"(cosh77 — 1) 



(cosh 770 — 1)2 
sinh 7/ 



(cosh 77 — 1)^ 
sinh 77 



cosh 77 — 1 



3q 



(2 cosh 7/ + 1) (sinh 77 — 7/) 
(cosh 77 — 1)2 



sinh 7/0 



9(w + 3a)c2/3 sinh 



Vo 



4(cosh77o - 1)2 



(cosh 77 — 1)2 (cosh 770 — 1)2 
4 sinh 77 9c^/2(2cosh77 + 1) 



9c3/2Q,(-2cosh7/ + 1) 
2(cosh77 — 1)2 



cosh 77 — 1 2(cosh77 — l)2r" 



(3.83) 
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where 



(C, - -C3) 



(3.84) 

(3.85) 
(3.86) 



a and modes 



and sinhr/ — rj = %c^/'^{r " — p)/2. 

The regularity condition at r = and r — oo impHes that modes with lo 
with < w < a are aUowed. In this case, this fixed point is a repeller up to the zero mode because 
ah other modes which satisfy the regularity condition have a positive w. Note that a suitable linear 
combination of the A^ and A J modes will have an asymptotic behavior k, 7.'^-2q Therefore the 
possible value of lo ranges from to 2a. The flow of RG in the vicinity of the fixed point is shown 
in Fig.d). 

(ii) case 2 (p = 0): 

The gauge condition which we use is, 



(3.87) 



By using the above condition Eq.(3.87), we obtain the physical perturbation 



{L) 

6p 



P a 



^(cosh?7-l) 

sinh rj 
cosh 77 — 1 



sinh// 

cosh — 1 



A 



2(coshry — 1) 



3 sinh 77(sinh rj — rj) 
cosh 77 — 1 



(2 cosh 77 + 1) (sinh 77 — rj) 



(cosh 77 — 1) 



(w + 3Q;)r" sinhr; 9q!C^/'^(2 coshr] + 1) 



cosh r] — 1 



2(cosh77 — 1) 



(3.88) 



where 



A+=ci/2c3. 



(3.89) 
(3.90) 



The regularity condition at r = and r = 00 implies that no modes are allowed. 

From the linear perturbation analysis around the fixed point, we see that the long-time behavior of the 
spherically symmetric dust universe is separated into two types. One is the case that the fixed point is 
a repeller. In this case, the Tolman-Bondi solution with self-similarity dose not play an important role 
in expanding universe because this fixed point is unstable and the spacetime will diverge from this fixed 
point. In the other case, the fixed point has both relevant and irrelevant modes. Although this fixed point 
is not a repeller, it has continuously many relevant modes. Thus it is not as straightforward as in the 
case of gravitational collapse Q to extract the long-time behavior of the universe, because it is sensitive 
to the initial condition and therefore we cannot uniquely predict the outcome. In the final section, we 
briefiy discuss how to treat the fixed point which has many relevant modes of the perturbation. 



IV. SUMMARY AND DISCUSSIONS 



We considered the spherically symmetri c but i nhom ogeneous universe filled with dust, where the Ein- 
stein equations have scale invariance Eqs.(3.39)-(3.42) and applied the renormalization group method to 
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study its long-time asymptotics. The fixed point of the RG transformation is a self-similar solution with 
scale invariance of the Einstein equations. In order to study the flow of the RG around this flxed point, 
the linear perturbation analysis is used. We impose the perturbation on the regularity at the boundary 
where the radial coordinate r equals zero or infinity. This boundary means that the area radius equals 
zero or infinity in the case of c = 0, on the other hand, in the case of c > 0, it equals a finite or infinity. 
The fixed point is the Tolman-Bondi solution with self-similarity, which includes the flat Friedmann uni- 
verse. The behavior of the fixed point is separated into two types. Both types have many relevant modes 
of the perturbation. 

The Tolman-Bondi solution with self-similarity is unstable against almost all spherical modes of linear 
perturbation. The spacetime will deviate from this fixed point. It is necessary to study the non-spherical 
mode of perturbation to say something more definite. In the cosmological problem, only the statistical 
quantities are meaningful if we think of comparison with observations. There are some works in the 
RG approach [p|-pO| on the universe which has a hierarchical structure. We may contemplate further 
development of the RG approach to cosmology formulated in the present work by introducing some 
kind of volume or statistical average for observables like energy density and the Hubble constant of the 
universe. The statistical concept is needed not only for comparison with observations but also for us to 
proceed further in the analysis of the RG equation because we have continuously many relevant (growing) 
modes around the fixed points. That is, the long-time behavior of the universe is sensitive to the initial 
configuration which we have no a priori control and we have to consider statistical likelihood of the initial 
values. 

We remark that the introduction of volume average in a finite region of the universe potentially intro- 
duces the scale invariance violation by hand because the exact scale invariance holds only for an infinite 
space. Note that in quantum field theories and statistical physics of the second order phase transition 
the scale invariance violations are hidden in the form of cut-off of the spectrum of physical modes. We 
shall elaborate our present observation in our future work. 

The self-similar solution given by Eqs.( 3.52 ) and ( 3.53 ) through the fixed point of the RG equation 
is essentially a function of t~^/"r — t^~~^l°'r jt^ which is roughly the fraction of physical distance to the 
hori zon scale of the Friedmann universe. Also note that in the case of non-linear diffusion equation 
Eq. ( ^.6[ ) implies the self-similar solution is a function of the ratio of the distance x to the diffusion length 
\Jt. In the both cases, the self-similar solution is a function of the distance in the unit of physically 
relevant time dependent scale. We believe this is a general phenomenon and the physical background of 
the RG equation which governs how dynamical variables deviate from the self-similar solution. 
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APPENDIX A: TOLMAN-BONDI SOLUTION WITH SELF-SIMILARITY 



In the spherically symmetric universe filled with dust, the most general solution of the Einstein equa- 
tions is the Tolman-Bondi solution 1 1 1 1 : 



A(r,i) 



(Al) 



B(r,t) 



2/3 



for Ci(r) = 



^ (coshr? - 1) , [t- Csir) = (sinhyy - r,)^ for Ci(r) > 

2^^(1-^03 7?), (^f - C3(r) = ^^W-^(r?- sin r;)j for C7i(r)< 



(A2) 
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^Pir.t)^^, (A3) 

where Ci(r), C2{r), and C^^r) are arbitrary functions of r and a prime denotes the derivative with respect 
to r. By taking Ci(r]_= c . Cyj r) — 4r"/9, and C3(r) — pr"", we obtain the Tolman-Bondi solution with 



self-similarity Eqs.(|3.59|)- (|3.67|) 



As for the calculation of linear perturbation, since we concentrate on the spherical modes of perturba- 
tion around a self-similar solutions it is enough to consider the linear perturbation of arbitrary functions, 
Ci(^), C2{r), and C3(r). The perturbed quantities, 5Ci{r), 5C2{r), and 6C3{r), can be expressed by a 
superposition of modes with different u and taken in the following form, 

(A4) 
(A5) 
(A6) 



By coordinate transformation of r. 



<5Ci(r) = 


cir^, 


SC2{r) = 








r > r 





(A7) 



where F{r) is an arbitrary function of r, we obtain the gauge mode of linear perturbation. This function 
F{r) also can be expressed by the superposition of modes with different lu in the form, 

F(r)=/r-+i. (A8) 
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